1 

Time as a Statistical Variable and Intrinsic Decoherence 

Rodolfo Bonifacio 
Dipartimento di Fisica, Universita degli Studi di Milano 
INFN and INFM, Sezione di Milano, Via Celoria, 16, 20133, Milano, Italy 



P.A.C.S.: 03.65 

We propose a novel approach to intrinsic decoherence without adding new assumptions to 
standard quantum mechanics. We generalize the LiouviUe equation just by requiring the dynamical 
semigroup property of time evolution and dropping the unitarity requirement. With no approxi- 
mations and specific statistical assumptions we find a generalized LiouviUe equation which depends 
on two characteristic time Ti and T2 and reduces to the usual equation in the limit Xi^Xi— >0. 
However, for Xi and T2 arbitrarily small but finite, our equation can be written as a finite difference 
equation which predicts state reduction to the diagonal form in the energy representation. The rate of 
decoherence becomes faster at the macroscopic limit as the energy scale of the system increases. In 
our approach the evolution time appears, a posteriori, as a statistical variable with a Poisson-F 
function probability distribution as if time evolution would take place randomly at average intervals 
T2 each evolution having a time width Xi. This view point is supported by the derivation of a 
generalized Tam Mandelstam inequality. The relation with previous work by Milburn, with laser and 
micromaser theory and many experimental testable examples are described. The agreement with 
recent experiments on damped Rabi oscillations is discussed. 
1. Introduction: 

The existence of coherent superposition of states is the basic reason for many paradoxical 
aspects of quantum mechanics. The evolution from a coherent superposition state to a statistical 
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mixture is called decoherence. This is a central problem for measurement theory and for the classical 
limit of quantum mechanics at the macroscopic level. The Schroedinger cat, which can be in a 
superposition of states dead or alive, or a macroscopic particle which can be in a superposition of 
"here" and "there" are typical example of paradoxes whose interpretation is stiU controversial. A 
superposition state gives non zero off-diagonal elements of the density operator, which give rise to 
quantum interference and non-classical correlation effects, von Neuman postulated the reduction to 
the diagonal form as a result of a measurement. However, this reduction remains mysterious because 
it cannot be described by a unitary Hamiltonian evolution and does not clarify how and when state 
reduction takes place and what the underlying dynamical process is. Essentially two approaches to 
decoherence have been proposed: the most widely used [1,2] takes quantum mechanics as it is, 
appealing to dissipation due to the interaction of the system with the environment or, equivalently, 
with the measuring apparatus, tracing at the end, over the many degree of freedom of the 
environment. The decoherence of the system is usually described by Master Equations (ME) which 
are derived using reasonable statistical assumptions and specific hut arbitrary models to simulate the 
environment and its interaction with the system. Furthermore, many approximations such as 
perturbative expansion, Markov approximation, etc. are used. From a conceptual view point, it 
appears rather peculiar that one has to invoke the environment, dissipation or the measuring 
apparatus to conclude that the cat must be dead or alive and a big particle must be here or there and 
not in a superposition of these possibilities with interference between them. It is like saying that the 
particle is going to be localized here as a result of is friction, and that the macroscopic limit is due to 
in the fact that the damping becomes stronger in this limit. 

A minority of scientists look for a modification of standard quantum mechanics to include 
decoherence as something intrinsic, i.e., not related to any specific model of interaction or 
entanglement with the universe and without perturbation and expansions and Markov 
approximations. The most widely known solution has been proposed by Ghirardi, Rimini and Weber 
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(GRW) [3]. They modify phenomenologically the Schroedinger equation by adding nonlinear terms 
to the Harmltonian so that the system, at Poisson distributed times, undergoes a sudden localization. 
Their model contains two unspecified parameters: the frequency and the spatial extent of 
localization. Other similar models have been proposed [4,5]. However, a common feature of the 
GRW model and of similar models is that the energy of the system is not preserved. Here, we look 
for an intrinsic mechanism for decoherence, which acts only on the coherence and not on the energy 
of the system. 

A solution in this sense has been proposed by the author in 83 [6] assuming a finite difference 
LiouviUe equation with a time step x (cronon) which has been related to the time-energy uncertainty 
relation. In ref. [7] it has been shown that this finite difference equation is equivalent to a semigroup 
master equation of the Lindblad form. 

Recently, MUburn [8] has proposed a modification of the LiouviUe equation assuming that: 
"the system does not evolve continuously under unitary time evolution, but rather on a stochastic 
sequence of identical unitary transformations" according to a Poisson distribution. However, in the 
Milbum theory time is a parameter, as in standard quantum mechanics, and the Poisson distribution 
is assumed as the limit of a Bernoulli distribution. 

In this paper we generalize the LiouviUe equation without any specific statistical assumption 
and without using any specific model for the interaction with the environment. We require only the 
semigroup property of the time evolution, dropping the unitarity requirement. We obtain an 
expression for the density operator which, a posteriori, can be interpreted as j/the evolution time is 
not a fixed parameter but a statistical variable whose distribution is that of the waiting time of "Une 
theory", i.e., a F-Poisson distribution. Our time distribution function contains two characteristic 
times which appear naturally in the theory as scaling times. They appears as the analog of time in the 
GRW localization of space. To be precise, time evolution appears as a series of random events in 
which T2 is the average time step between two evolutions and Xi is the time duration of each 
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evolution. It is the same as saying that T2 is the average interval between the arrival of two "clients" 
at the counter in a bank and Xi is the time each "client" spends at the teUer. Therefore, Xi and X2 have 
very different physical meaning and in general Xi < X2. Our time evolution law, for Xi = X2 ^ 
recovers the Liouville equation, whereas for Xi and X2 arbitrarily small but finite, results in a 
irreversible state reduction to the diagonal form in the energy representation. This state reduction has 
been postulated in ref.[9]. Here it is dynamically derived. Assuming x = Xi = X2 sufficiently small, we 
obtain the Milburn ME [8]. We suggest an application of our formalism to continuous measurement 
theory: X2'' is the observation rate and Xi is the time duration of each observation. Unlike the 
treatment of ref. [3], ours has energy as a constant of motion. Our evolution equation can be written 
in the form of a finite difference equation with time step X2. Therefore, according to our formalism, 
one cannot give a continuos and instantaneous description of time evolution but only a description by 
"quantum jumps" with time steps given by the cronon X2. From the finite difference equation we 
obtain a generalized Tarn Mandelstam inequality which connects Xi and X2 with the time-energy 
uncertainty relation. Extending our formalism to a non Hamiltonian system we derive a 
generalization of the ScuUy-Lamb ME for laser and micromaser. Here Xi and X2 have a precise 
physical meaning because the micromaser model has been extensively investigated theoretically and 
experimentally to measure the quantum decoherence of an e.m. field prepared in a Schroedinger cat- 
like state, as one approaches the macroscopic limit (increasing the photon number) [10 - 13]. We 
apply our formalism to many examples which can be experimentally tested. In particular we derive an 
extra diffusion term in the position spread of a free particle; we describe localization of a free particle 
prepared in a Schroedinger cat-like state of two different positions. We demonstrate the existence of 
an intrinsic linewidth for a single mode e.m. field. Finally, we describe decoherence for a spin 
superposition state in a magnetic field and cancellation of EPR correlation under conditions in which 
they should not disappear according to standard quantum mechanics. Furthermore, intrinsic damping 
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of Rabi oscillations in a two-level system, in agreement with recent experimental results [10, 14] is 
discussed. 

2. The finite difference Liouville equation: 

The unitary time evolution of a quantum system is generally described by the Liouville- von 
Neumann equation 

(1) -^ = -iLp(t), 

where p is the density operator, L is the Liouvillian super-operator Lp = — [h,p1 , and H is the 

n 

Hamiltonian. The formal solution of (1) can be written as 

i i 

(2) p(t) = e-^'Po=e"^"poe^"'- 

In (1) and (2) time t appears as a parameter, not an "observable" or a statistical variable as, for 
example, are position, momentum and H. Equation (1) and (2) can be expressed in the energy basis 
Hln>=Eln> as: 

(3) P„,m = -i«n,mP„,m 

(4) P„,.(t) = p„,JO)e-'"'»»' 

where (0^^=(Ej^-Ej^)//z. The degenerate case can be included in this notation by assuming that 
the same states ln> belongs to the same eigenstate En. In the following treatment, in case of 
ambiguity with super-operators, one can refer to matrix elements replacing L with a)t,m as it appears 
fromEqs.(2) and (4). 

Let us deiine a generalized density operator p defined as 



(5) p(t) = jdt'P(t,t')p(t'); t,t'>0 



where P(t,t') is a function still to be determined. In particular if P(t,t') = 5(t-t'), p(t) = p(t) , so that 
p(t)is a generalization of p(t) , if P is unspecified. Using Eq.(2) we can write 
(6) p(t) = V(L,t)p(0) 

where 



(7) V(L,t) = jdt'P(t,t')e" 



iLt' 



We now determine P(t, t') imposing the following conditions: 

i) p(t) = p"(t)>0; Trp(t) = l 

ii) V(t + t') = V(t)V(t') 

Condition (i) identifies p(t) as a density operator. Condition (ii) is the so called semigroup property 

which ensures translational invariance of the initial condition i.e., p(t + 1') = V(t + t')p(O) = 

V(t)p(t') = V(t')p(t) . Note that (i) and (ii) are satisfied by the usual LiouviUe operator V(t) = e"'^' . 

We just drop the request of Unitarity = V"^ Condition (1), using (5) and taking the trace on 
both sides implies that 

(8) jdt'P(t,t') = l; P(t,t')>0 



where we used Trp / 1') = 1 . Therefore, Equation (8) defines P(t,t') as a probability distribution 
function; P(t,t')dt' can be read as the probability that the random variable t takes a value between t' 
and t'+dt'. Therefore, p(t) and V(t) appear as an average value respectively of p(t) and V(t) . The 
semigroup property can be generally satisfied assuming 

(9) V(L,t) = [Vi(L)r^ 
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where T2 is a scaling time. Note that also the Liouville operator V(t) = e can be written in the 

iLt 

form (9) taking Vj = e~ ^ (unitary). We now determine the most general form of Vj(L) compatible 
with the previous requirements. Let us use the F function integral identity, 

(10) (A+iB)-'' = jdX^e-^'e-^' ; 

which is valid for A>0 and k>0. Let us identify the following: 

Vi~'=A+iB; k = t/T2; 'k = \:ix, 

where Ti is a scaling time generally different from T2. In this way by imposing the consistency of 
Eq.(lO) with (7) for aU L and the normalization condition (8), we obtain B = LTj and A = 1 so that 

(11) p(t) = v(t)p(o) = - — ^^p(o) 

(l + iLxJ 



(12) p(t) = -— ln(l + iLx^)p(t) 

and P(t, t') is the r distribution function: 



(13) P(t,t')= ' 



1 e"'''^' f t' 



(t,t'>0) 



1 r(t/xj 

This argument can be supported by considering P(t,t')=0 for t'<0 and using the unicity of the Fourier 
transform. Equations (11), (12) and (13) are the basic result of our paper. In particular, equation 
(11) provides a generalized form of the density operator and of the time evolution law which has 
been obtained just by requiring the density operator properties (i), the semigroup property (ii) and 
dropping the unitary requirement. Therefore, we are not adding new elements to the formalism of 
quantum mechanics but, on the contrary, we are reducing the basic assumptions. Taking 
Tj = T2 — > in Eq.(l 1) one obtains the Liouville limit and P(t,t') = 6(t - 1') . However, as will be 
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shown later, for Ti and T2 arbitrarily small but iinite, as we shall see, one irreversibly approaches the 
diagonal form in the energy representation. 
3. The Finite Difference Equation: 

Equation (1 1) in the Liouville limit Tj = ^ recovers again Eq.(l). When Xi and T2 are iinite, the 
second order expansion of Eq.(l 1) gives 



where L^p = [h,[h,pJJ . Equation (14) with Xi = T2 gives the well known "phase-destroying" Master 

Equation (ME), deduced by many authors using a reservoir interaction model [2] or specific 
statistical assumptions [8] and used in Quantum Non Demolition (QND) measurement theory [11]. It 
is straightforward to show fromEq.(ll) that p(t) obeys the following finite difference equation 



where H = Hx^ . Equation (15), for = — > , again gives the Liouville equation, whereas for 
Xj = X2 = X it reduces to an equation proposed a long time ago [6] to describe irreversible state 
reduction to the diagonal form. The difference here is that now it has been derived under very 
general assumptions and two characteristic times appear, with very different physical meanings. In 
ref [7], it has been shown that for Xj = Xj = x , Eq.(15) is equivalent to a semigroup ME of the 

Lindblad form. The same considerations apply here to Eq.(15), taking X2 = x and substituting H with 
H . In this way Eq.(15) is formally identical to that considered in ref [6]. 

Equation (15) shows a very important feature. If p(t)is a solution of (15), then f (t)p(t) is 
also a solution, provided f (t + Xj) = f(t) . Therefore, p(t)is uniquely determined only within the 
"cronon" Xi, i.e. for time intervals t = kXj , with k integer. This fact implies a redefinition of p(0) 
which in the standard description is the density operator determined at some instant t=0. This is 



(14) 




(15) 
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clearly an artifact, because the possibility of an instantaneous measurement of a complete set of 
observables to determine p(0) at the instant t=0 appears as a mathematical abstraction. Our iinite 
interval description of time evolution, which follows from Eq.(15), appears much more realistic 
because p(0) can be interpreted as the density operator determined in a cronon interval T2. The 
evolution of Eq.(15) on the "time grid" gives p at later time intervals, k=t/X2 > 1, and can be 
parametrized as 

(16) p(k)-p(k-l)=iLXiP(k), k>l 

or equivalently p(k + 1) = (l + iLx^ ) ' p(k) ~ e p(k) , for Xi smaU enough. It is as j/time evolution 

occurs discontinuously in "quantum jumps" [12] spaced by intervals X2 with each evolution occurring 
within a width Xi. Each "jump" is given by a unitary time evolution only if Xi is small enough. 
Accordingly equations (7) and (11), fort/X2=k>l, with r(k) = (k - 1) ! can be written as: 

(17) V(k) = |dt'P(k,t')e"^'' ; 



and 

''*''''=^:= 

where P(k, t') can be interpreted in two ways: i) as the weU known Poisson distribution in k, or ii) as 
the F-distribution function in the continuous variable t'. Unlike previous treatments [8, 13] we adopt 
the latter. Equations (16) and (17) can be interpreted in terms of the waiting time statistics for k 
independent events. According to (17), time evolution is made up of random "events", i.e., unitary 
time "evolution". The probability density for k=t/X2 events to take place by a time t' is given by 
Eq.(18). In particular, X2 is the average interval between two "events" ( Xj^ is the rate) and Xi is the 
time width of each event. We propose an interpretation of our formalism in terms of a continuos 
measurement theory: x'^ is the observation rate, Xi is the time width of each observation. As we 
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shall see, this interpretation has a precise meaning in cavity QED measurements as well as in 
micromaser and laser theory [13]. In our formalism the interaction with the measurement apparatus 
is described by two characteristic time Ti and T2: this characterization is intrinsic and model 
independent in the sense that it does not depend on the way the measurement is carried out or on the 
detail on the measurement apparatus. Because the Harrdltonian is a constant of motion, our 
formalism can be applied, though not exclusively, to a QND measurement. According to this 
interpretation one obtains a dynamical definition of state reduction in agreement with the von 
Neumann state reduction. In fact, in general, in the energy representation Eq.(ll) becomes 

(19) P„,.(t) = 7 ^^P„,.(0) = . . ^ = e-^^-'e-^^»-'p„, JO) 

where 

(20) Y„,„ = (l/2Tjln(l + <^T^); v„,„ = (l/x Jarctga)„,„Ti . 

Note that irreversibility is obtained also in the limit Xj — > and — > , provided that xf /Xj is 
finite, which describes an irreversible evolution. In general, the decoherence rates Yn,m increases as X2 
becomes smaller or Xi becomes larger. Assuming, for simplicity, non degeneracy, for En =Em , 
Therefore, tOh,m=0 and p„ „ = Pn „(0) , so that the energy is a constant of motion, whereas for n^m. 



Pn,, 



with a rate Yn,m. Therefore, p(t) — > ^p„„(0)|n)(n| i.e., p approaches the stationary 



diagonal form. Then a pure state remains a pure state if and only if it is a stationary state; otherwise 
the system wiU evolve to a statistical mixture. Note that, at the microscopic limit, when the energy 
scale becomes very large, the decoherence time l/y becomes very short. The rate constants Yn,m and 
the frequencies Vn,m have different behavior depending on whether co^^Xj » 1 or whether 
to„^Xi « 1. If co„^Xi « 1 one has 
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(21) 




and 




to 



where we have omitted reference to m and n. This is also the result one obtains by the phase 
diffusing ME (14). Therefore, this ME is valid only if CO^^Xi « 1 , for aU n and m. This assumption 

can be made only if the spectrum is bounded so that it cannot be applied to a simple harmonic 
oscillator, as is usually assumed [11]. In the opposite limit, to^^Xj » 1, one has 



Therefore, for large to there is a slow logarithm dependence of the decoherence rate on the energy 
separation whereas v becomes independent of to. We are now in a position to interpret the behavior 
of the frequencies tOi,m as given by Eq.(21). If tOi,m Xi « 1 the decoherence time tD=l/Y is much 
larger than Xi. Therefore, it takes many Xi time step for decoherence to occur. Regarding the 
frequency, we note that using Eq.(21), e'^' = e'"" , where t = (t/Xj )xj is the "effective evolution 
time" given by the number of evolution steps t/X2, times the width of each evolution step, Xi. If tOi,mXi 
» 1 one has, from Eq.(22), a decoherence time to = l/y - '^2 /ln(caXi)«X2. This implies "immediate 
decoherence", i.e., von Newmann state reduction after the first evolution interval X2. In this case of 
strong decoherence, the "frequency" loses its significance as it is smaller than the linewidth. We 
emphasize that the basic point of our treatment is that in our treatment time appears as a statistical 
variable with a distribution function P(t,t'), given by Eq.(13) and (18), which for k=t/X2=l, is a 
simple exponential. However, for t/X2 » 1, P(t,t') is a strongly peaked function with a mean value, 
<t'>, and dispersion, o, given by 



(23) 



<t'>=t = (t/x2)x 



1' 



(24) 
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Note that <t'> does not coincide with t but with the effective "evolution time" t . The dispersion a 
scales as Vt , like in a diffusion process. According to the previous interpretation, the dispersion of a 
appears as the dispersion due to k=t/X2 statistically independent "events" times the width of each 
event Ti. For k = 1, o assumes its minimum value Ti. Therefore, Ti appears as an "inner time" [15], 
i.e., the intrinsic minimum uncertainty of the time evolution. The relative dispersion, a/ < t'>= 
(t / X2 ) , goes to zero as the number of time steps, t/Xi, goes to infinity. Furthermore, for t/Xi »1, 

P(t, t') can be approximated by a Gaussian in t' with mean value given by Eq.(23) and dispersion 
given by Eq.(24). 

4. Comparison with previous work 

The difference between our theory and previous descriptions of intrinsic decoherence [8] can 
be summarized as foUows: In ref. [8] it is assumed that i) the system evolves under a random 
sequences of identical unitary transformations, and ii) the probability of n transformations in a time t 
is given by a Poisson distribution so that (using our notation) 

(25) p(t) = £p(n,t)e-""^'p(0) 

(t/xj° 

where p(n,t) = j — e ' . In ref. [8] Xi = X2 = Y is interpreted as a "fundamental time of the 



universe". Summing up the series (25) one obtains p(t) = exp 



— (e"'^'' -l)\(0) so that 
[^2 ] 



(26) ^(t) = — (e-^^' - l)p(t) = i-(e-^^''''p(t)e"'''''' - p(t)) 



'C2 '^i 



Equation (26), with Xi = X2 = is the basic result of ref. [8]. If one applies Eq.(26) with this 
assumption to calculate the average value of the annihilation operator of an harmonic oscillator, one 
obtains easily < a >= Y(e~'"'^ - 1)< a > . From this equation Milburn [8] infers that there are particular 
frequencies (x)=2mvY where the oscillator is "frozen", i.e., < a >= . As we shall see this unusual 
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behavior does not occur in our treatment because we always have damped amplitude for all 
frequencies. We note that the same method has been adopted in Laser-Micromaser and cavity QED 
theories [13], where the formal substitution: 

(27) e"^''=>M(Ti) ; L=>(i/TjlnM 

have been made in Eq.(25). Here M is a non unitary operator well known in laser-micromaser theory 
described in ref.[13]. In this way Eq.(26) becomes the well known Scully-Lamb ME 

(28) P = (1/tJm(Ti)-1)p 

where = r is the rate of injection of the atoms and Xi is the atomic interaction time. The basic 
differences between our approach and Milburn approach are: 

i) we do not make any statistical assumption: randomness in time evolution appears naturally as an 
interpretation of our results, i.e., of Eq.(7) and (13). In particular the Poisson-F distribution, 
Eq.(13), is not assumed but derived. 

ii) the choice Xi = X2 implies a severe restriction on the statistical interpretation of the theory. In fact, 
in our treatment, as well as in the Scully-Lamb theory, Xi and X2 have very different physical 
meanings. 

ui) Equation (25) assumes that the number of evolution transformations in a time t is a random 
variable and t is just a parameter as in the Schroedinger equation or in the Liouville equation. In 
our theory the waiting time t for n evolutions to occur is the random variable. 

The two view points appear similar but are basically different. In fact, they lead to different results, 

as one can see by comparing Eq.(26) and Eq.(12). The first can be derived from the second as 

follows. Using the integral identity: 

-X 

(29) ln(l + ix) = f d^i V (1 - e"'^' ) 

^ 

and taking x = Lxi and A,=t' /Xi, Eq.(12) can be written in the integral form 

13 
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(30) 



mt) = ^jdt' 

'^2 







"e-'^''-f 












t' 


J 





p(t) 



This equation is completely equivalent to Eq.(12). It can be shown that if the term in square brackets 
is slowly varying on a time scale Xi, it can be taken out of the integral with t'= Xi. Because the term 
in the round brackets is normalized, one obtains the Milbum ME (26). Therefore, Eq.(26) can be 
obtained as an approximation of Eq.( 12), if Xi = X2 = is small enough. In the case of an harmonic 
oscillator this implies the assumption co/y «1, which is inconsistent with the condition for freezing 
03/'Y=2oK. Note that by making the substitution (27) in Eq.(12) and (30) one obtains a generalized 
laser-micromaser ME: 



(31) 



P(t) 



jdt'-e-''^' 



'^2 '^1 



M''^^--l 
t' 



p(t) = -— ln(l-lnM)p(t). 



Again if Xi is small enough, Eq.(31) reduces to the Scully-Lamb ME (28). Therefore, Eq.(31), with 
no restrictions on Xi, appears as a generalization of Eq.(28). As already mentioned, in the context of 
cavity QED [13] the two characteristic times have a very precise physical meaning. It is interesting to 

note that the finite difference equation (15), with the substitution (27), reads — — = 



-InM 



p(t) . Replacing the RHS with the derivative, one obtains 



(32) 



dp 
dt 



-InM 

V ''2 J 



which is the equation of a micromaser with regular injection [13]. Finally, dissipation can be taken 
into account by adding the proper dissipative term of the Lindblad form to Eq.(12) or (31). 
The possibility of applying our formalism to cavity QED is of particular relevance to the quantum 
theory of measurement. Cavity QED has been extensively investigated both theoretically and 
experimentally [10,13]. 
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5. Generalized Tarn Mandelstam Relation: 

From the finite difference equation (15), and defining the mean value of an observable A in the usual 

way, A = Tr p A , we obtain 

(33) A(t)-A(t-T,) _|^ 

Let us apply Eq.(33) to the one dimensional motion of a particle with Harmltonian 

H = p^ / (2m) + V(x) . Taking A=x and A=px we have (x(t) - x(t - X2 ) / = p^ (t) and 

(Px (t) ~ Px (t - '^2 ) / '^1 = F(x) > where F(x) = -dV / dx . This is the finite difference version of the 
Herenfest theorem. The usual form is regained in the continuos limit = ^ , whereas the 



classical limit is F(x) = F(x) . Note that the two limits are independent. Therefore, taking only the 
classical limit one obtains finite difference Hamiltonian equations for x and p^ . 

We now derive a generalized Tam Mandelstam relation. Using the general uncertainty relation for 
A and H we can write 

;z|aa| 



(34) a(A)o(H)>i|[Ajr| 



2 Xj 

where AA = A(t) - A(t - Xj ) and we have used Eq.(33). Equation (34) can be written in the form 

(35) x^a(H)>;//2 ; x^^^^. 

|AA|/Xi 

This equation appears as a generalized TM inequality and reduces to the usual form in the limit Xi=X2 



0. In fact in this way one obtains x^ = a(A) / 
inequality 



A 



However, going back to Eq. (34) one obtains the 



AA 



X, 



(36) 

a(A) X 
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where TE=^/2a(H) is the intrinsic inner time [15] of the system. Because inequality (36) is valid for 
any observable, we can conclude that if Ti < Te then lA A l< a(A) necessarily follows. Therefore, no 
appreciable variation occurs for any observable in the cronon time Xi. We can say that if Xi< Xe one 
has a quasi continuous evolution of the system even using the discrete time description. On the other 
hand, Eq. (36) states that if X2 is such that lA A l> a(A), it follows that one must have Xi > Xe. 
Therefore, the ratio Xi / Xe rules the rate of change of the state of the system within the time interval 
X2 between two evolutions. The choice Xi = Xe = fi /2a(H) corresponds to the maximum possible 
value of Xi which guarantees the quasi continuous evolution which is commonly observed. This 
choice is clearly meaningless in case of a non Hamiltonian system, as in Eq.(31). Furthermore, from 
X2 > Xi = Xe, one would obtain X20(H) >h/2, giving a precise and intrinsic meaning to the time- 
energy uncertainty relation. 
6. Calculation of Physical Quantities: 

Equation (5) allows one to calculate aU physical quantities by just performing a time integral of the 
usual expression. Let us give some examples: The mean value A(t) = Trp(t)A can be easily 
obtained using Eq.(5) as 

(37) A(t) = jdt'P(t,t')<A>, 



where 

(38) <A>,= Trp(t)A 

is the usual expectation value. Therefore, a constant of motion remains a constant of motion, 
whereas oscillating quantities, as e'"' , are damped with a rate constant y = (l/ 2X2 )ln(l + co^xf ) 
similarly to pn,m as given by Eq.(19). Furthermore, from Eq.(5), 

(39) (x|p(t)|x'} = jdt'P(t,t')(x|p(t')|x') . 
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The same relation holds for matrix elements in any representation. However, Eq.(39) is of particular 
relevance. In fact, for x = x' it gives the position probability density, P(x,t) . In particular, if the 

initial state is a pure state, | , then (x |p(t')| x) = |v|/(x, t')|^ . Furthermore, because the 

Wigner function is related to the Fourier transform of ^x|p(t)|x'^ it can be obtained using the same 

integral relation as in Eq.(39). It can be easily shown that the same relation is valid for the Glauber 
P-function. We now make some testable application of our formalism. 
7. Harmonic oscillator: 

In this case the photon number is a constant of motion. On the contrary the amplitude <a> goes to 

zero. In fact, because < a >j=< a >o e"'"" , from (37) and (8), one obtains a(t) = ^- — r-r^ < a >q= 

(l + icoii ) ' 

< a >Q e"^e"'^' , where y = ^^^^^(l + ^^'^i ) ^ ~ "arctgcoXj . Therefore, |a(t)| goes to zero as t 

goes to infinity. The behavior of y and v is formally the same as that of yn,m and Vn,m described in 
Eq.(20). Therefore, the same discussions apply just abolishing the index n,m. Moreover, in general, if 
one assumes that the initial state is a coherent state, lao>, the state at time t will be the coherent state 
I a(t)} with a(t) = age"'"" . Therefore, 

(40) p(t) = ^e-' e'"^°-"'1n)(m| 

Vn!m! 

where n = |ao| . The oscillatory behavior of the diagonal elements implies that as t — > oo , 

p(t) ^e~"(n)" /n!. This happens because p„,„(t) = e"^°"' with y„^ = (l/T2)ln[l + a)^(n - m)^Tf J, 

n 

where Eq.(21) has been used. Note that the decoherence rates increase as one increases the distance 
between energy levels. For n and m such that a)|n - m|Tj « 1 , expanding the logarithm to the first 

order, y„ ^ oc n - m , as has been already derived by other authors [11] with the ME (20) whereas if 
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In-ml is large enough Yn,m has only a logarithmic dependence on In-ml and the ME (14) does not apply. 
The diagonal density operator (40) has been introduced by Glauber [16] to describe the steady-state of 
a laser well above threshold using a random phase assumption for ao. Here, the approach to the 
diagonal form is dynamically obtained even if the phase of ao is perfectly determined. Note that the 
amplitude damping should be experimentally observable by proper homodyne detection or, more simply, 
by measuring the linewidth. This intrinsic linewidth y should be observable only if it dominates on all 
other linewidths, i.e., if Ti is large enough and T2 is smaU enough. In the opposite case one obtains a 
lower limit for X2 and an upper limit for Xi. 

8. Free particle spread: 

Let us assume to have a free particle in a minimum uncertainty packet, (x|p|x^ = |\|/(x,t)| = 

(l/V27Uo(t))e"^''"^''^''''''''*'\ where <x>^=Xo + vt, a'(t) = +o't' with v = p/mand 

= Op / m . Using the integral (39) with x=x' and a Gaussian approximation for P(t,t') one obtains 
again a Gaussian for P(x,t) centered on x(t) = Xq + vt and spread 

(41) ^<5l+(5l~e+wh,t 

where t is the reduced time t = tXj / . If Xi = T2 , this expression coincides with the one derived in 

ref. [8]. Assuming v » Cy, this term is dominant, if t<(v/a^) Xj. So, if v/Oy and X2 are large 

enough, the diffusion term should be observable. 

9. Space Localization in a Superposition State: 

Let us consider a free particle prepared in a superposition state of two minimum uncertainty wave 
packets, centered at different positions xi, and X2. The wave function at time t is given by 

(42) \|/(x, t) = (1 / V2)(v|/i (x, t) + V|/2 (x, t)) 
where 
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(43) ¥j(x,t)= ,^^ ,.2 """'"''■^^^'"^ 



V ^^x ) 



j=l,2 



Here, for simplicity, we have neglected the Schroedinger spread, by assuming to be small enough. 
This is certainly the case at the macroscopic limit, because =h/2ma^ , where m is the mass of 
the particle. From Eq.(42) we have: 

(44) P(x,t) = |vi/(x,t)|' = (l/2^\|/^(x)|' + (l/2^\|/2(x)|'+Int 
where Int = |\|/i(x)||\|/2(x)|cosa)t, with 

(45) (0 = c^xD/2cl ; 

with the assumption that Xi = -X2 = D/2. Therefore, in standard quantum mechanics the interference 
term is oscillating in time. In contrast, in our formalism the interference disappears because, using 
Eq.(37) and (8), the average of cos(tot) is exponentially damped at a rate given by 

(46) Y = (l/2T2)ln(l + a)'Ti), 

where CO is given by Eq.(45). The interference disappears in a decoherence time to = Y Note that if 
cori«l one can expand the logarithm to the &st order, obtaining a decoherence rate proportional to 
of, i.e., to the square of the distance between the centers of the two packets, in agreement with 
previous treatments. The interaction with the environment and/or with the measuring apparatus 
appears quite naturally in our formalism via the two characteristic times without specifying any 
interaction model. If we use the uncertainty principle to eliminate , Eq.(45) can be written as 

(47) (xi = l6mE^Dx/h\ 

where E = (l/2)m < v^ > . This expression, together with Eq.(45), clearly shows that at the 
macroscopic limit, when m and E are very large, or at the classical limit — > , the decoherence rate 
becomes extremely fast. Let us remember the fundamental and mysterious meaning of the 
interference term. Its presence implies that to find the particle around Xi or around X2 are not 
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mutually exclusive events. Otherwise, this probability would be, by definition, the sum of the 
probabilities. Therefore, the interference term, independently of its magnitude, says, roughly 
speaking that the particle is in a "schizophrenic" superposition of here and there. This is substantially 
the Schroedinger cat paradox, if xi and X2 are replaced by "dead" or "alive". Quoting Fe5Timan, 
"interference is the only mysteries of quantum mechanics" [17]. We call decoherence or state 
reduction, the suppression of the interference term so that the particle can be found at Xi or X2 with 
50% classical probability, but without being able to predict where. We think that this is the "best" 
localization quantum mechanics can describe without strongly changing its formalism. Other authors 
[3] introduce non linear terms in the Schroedinger equation to have not only the suppression of the 
interference term but also the reduction to a single wave packet which, in the Copenhagen 
interpretation should take place only when one really observes the particle. Accordingly, only the 
suppression of the interference term should be intrinsic. Many authors [1,18] describe decoherence 
as due to dissipation, which is like saying that the particle is here or there because it is going to be 
stopped by friction or the cat cannot be dead and aUve because it is going to dissipate its energy in 
the universe and die. It is clear that dissipation is sujficient for decoherence. But it is not clear that it 
is necessary for decoherence. The interference does not disappear at x=0 because there the 
interference term is stationary. For this reason our decoherence mechanism is not present in the case 
of the stationary two slits interference experiment between particles having the same energy or 
momentum. In our case the time dependence is due to the finite velocity spread. The decay rate, even 
if it is zero exactly at x=0, is different from zero in a region whose extent goes to zero at the 
macroscopic limit. Since the probability is the integral of the probability density, the interference term 
goes to zero in a finite region around x=0. 
10. Decoherence in two level systems: 

Let us assume a spin system of neutral particles injected one at a time with velocity v, in a constant 
magnetic field Bo of extent L in the z direction. The particles are prepared in a eigenstate of Ox, 
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I ) = (I +) + 1 -)) . Therefore, \\\f)^=-j= (e'""' | +) + e"'""' | ->) . Here cOq = iiBq / is the Larmor 
frequency and |X is the Bohr magneton. The density operator at time t reads p(t) = | ^I')^^ | = 
^(pn +e'"''''|+)(-| + h.c), where p^ =(l/2)[|+>(+| + |->(-|] is the diagonal part. Note that for 

a)t)t=2n7t one has again the initial pure state. However, using our formalism the off-diagonal elements 
will be damped as e"^ where 

(48) Y = (l/2x2)ln(l + cOoXi) 

and t = L/v is the transit time in the magnetic field. Let us now insert, after the region L a Stem- 
Gerlach apparatus in the x direction and let us take tOtjL/v = 2aK with n arbitrarily large. According 
to standard quantum mechanics aU particles would go only in one direction. However, if yL/v is large 
enough, our decoherence will take place, i.e., 50% of the particle will go both directions. The non 
observation of this decoherence would put a upper limit for Xi and a lower limit for Xi- 
11. Intrinsic Damping of Rabi Oscillations: 

The same considerations can be applied to a system of two level atoms injected one at the 
time in a high Q resonant cavity and prepared in a Rydberg state so that the atomic and cavity decay 
time and the intrinsic decoherence are very long. This is the experimental situation in ref. [10]. If the 
field is in a n-Fock state and the atoms are injected in the excited state, the atoms wiU oscillate 
between the upper and lower state so that the population difference d will oscillate as d = cosHt , 

where Q. = g-Jn + l is the Rabi frequency and g is the one photon Rabi frequency accordingly with 
the Jaynes-Cummings atom-field interaction Hanultonian model. In our formalism these oscillations 
win be damped with a rate constant, y = (l/2x2)ln(l + ^2^Xj ). Therefore, if yt = yL/v»l, the 

population difference will approach the steady state value d=0 so that the population of the upper 
level will approach the value 0.5. This behavior has been indeed observed in ref. [10] also for the 
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vacuum state where dissipation losses are ineffective and in ref. [14] for a n-Fock state, in a different 
experimental situation. In particular, in ref. [14], it has been observed an increase of the damping rate 
Y with n, in "qualitative" agreement with (n+1)'^^. Let us point out that our damping rate has a 
logarithmic dependence on (n+1), which, if Xi is small enough, goes like (n+1). Whether this 
damping is due to experimental problems, as suggested by the authors, or can be attributed to our 
intrinsic decoherence mechanism deserves further experimental and theoretical investigation. 
Obviously the same intrinsic damping can be predicted for the Rabi oscillations in a NMR 
configuration. 

12. Cancellation of EPR correlation: 

Let us assume two spin- 1/2 neutral particles, say 1 and 2, traveling in opposite directions with 
velocity v and prepared in the singlet entangled state | V o ) = (|+i'~2)~|~i'+2)) with a constant 
magnetic field Bo in the z direction acting on the path of particle 1 for a length L. We have 
|\|/)^ = --^(e""°"^|+,-)-e~""°"^|-,+)), where coo is the Larmor frequency. The associated density 
operator reads 

/ I 

(49) p(t) = |x|/(t))(\|/(t)|= PD--e'"»'|+,-)(-,+ |-h.c 



where is the diagonal part. Note that for tOoL/v = 2n7U one has again the initial singlet state. 

Therefore, the presence of a e.m. field would become ineffective regarding EPR correlation under 
these conditions. However, in our formalism the oscillating of the off-diagonal terms disappears 
exponentially with y given by Eq.(48). Therefore, if yL/v »1 the EPR correlation should disappear 
even if cooL/v = 2n7t. Any evidence of correlations smaller than expected in the conditions described 
above, would be a strong support of our theory. 
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Conclusion: 

Using Anton Zeilinger [19] words: "When investigating various interpretations of quantum 
mechanics one notices that each interpretation contains an element which escapes a complete and full 
description. This element is always associated with the stochasticity of the individual event ... It is 
suggested that the objective randomness of the individual quantum event is a necessity of a 
description of the word. ... Yet it is also highly recommended to follow the guidance of the 
Copenhagen interpretation, that is, not to make any unnecessary assumption." 

In this spirit we have generalized the Liouville equation without introducing new assumptions 
to quantum mechanics. On the contrary, we have reduced the basic axioms of quantum mechanics 
dropping the unitarity condition and maintaining only the semigroup property of the time evolution 
operator. The new equation describes intrinsic decoherence giving irreversible state reduction to the 
diagonal form in the energy basis so it can be used to describe quantum non demolition 
measurements. In our treatment time evolution appears, a posteriori as a statistical variable obeying 
"the waiting time statistics" whereas in standard quantum mechanics as weU as in previous 
descriptions of intrinsic decoherence time is just a parameter. Our evolution equation can be written 
as a finite difference equation, which looks as j/time evolution is a random process with two 
characteristic times Xi and T2. The second time is the average interval between two evolution and Xi 
is the time width of each evolution. For Xi = T2 = x one recovers a previously proposed finite 
difference equation [6], which for x smaU enough and in a particular Umit, gives the ME for intrinsic 
decoherence proposed in ref. [8]. A generalized Tam Mandelstam inequality has been derived, which 
suggests a relation between Xi and the "inner time" of the system, /z/2a(H) [15]. This relation 

guarantees a quasi continuous evolution even with a finite difference equation. We have shown that 
our equation implies localization of a free particle prepared in a Schroedinger cat-like superposition 
of states with different position as the macroscopic limit is approached. The relation of our theory 
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with micromaser theory and experiments have also been discussed. Some testable examples have 
been described. In particular, we predict the damping of Rabi oscillations, when aU dissipation 
mechanisms appears to be ineffective, in agreement with the experimental result of Ref. [10] and 
[14]. 
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